Global modelling of X-ray spectra produced in O-type star winds 
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ABSTRACT 

High-resolution X-ray spectra of O-type stars revealed less wind absorption than expected from 
smooth winds with conventional mass-loss rates. Various solutions have been proposed, including porous 
winds, optically thick clumps or an overall reduction of the mass-loss rates. The latter has a strong impact 
on the evolution of the star. Our final goal is to analyse high resolution X-ray spectra of O-type stars with 
a multi temperature plasma model in order to determine crucial stellar and wind parameters such as the 
mass loss rate, the CNO abundances and the X-ray temperature plasma distribution in the wind. In this 
context we are developing a modelling tool to calculate synthetic X-ray spectra. We present, here, the 
main ingredients and physics necessary for a such work. Our code uses the most recent version of the 
AtomDB emissivities to compute the intrinsic emissivity of the hot plasma as well as the CMFGEN model 
atmosphere code to evaluate the opacity of the cool wind. Following the comparison between two for- 
malisms of stellar wind fragmentation, we introduce, for the first time in X-rays, the effects of a tenuous 
inter-clump medium. We then explore the quantitative impact of different model parameters on the X-ray 
spectra such as the position in the wind of the X-ray emitting plasma. For the first time, we also show 
that the two formalisms of stellar wind fragmentation yield different results, although the differences for 
individual lines are small and can probably not be tested with the current generation of X-ray telescopes. 
As an illustration of our method, we compare various synthetic line profiles to the observed O VIII A 18.97 
A line in the spectrum of ^Puppis. We illustrate how different combinations of parameters can actually 
lead to the same morphology of a single line, underlining the need to analyse the whole spectrum in a 
consistent way when attempting to constrain the parameters of the wind. 

Subject headings: stars: early-type - stars: mass-loss - X-rays: stars 

Introduction teractions with their environment. In fact, mass-loss 

plays a f undamental rol e in the evolution of massive 



The study of stellar winds is crucial for our under- stars (e g |Hirschi |[2008h as well as in the mechanical 



standing of the evolution of massive stars and their in- and chem i ca l feedback of these stars into the interstel- 



lar medium (e.g. iDopita Il2008f) . Over the last decade, 
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as the importance of wind fragmentation became clear. 

Direct observational evidence for the existence of 
clumps in the winds of O-type stars was obtained 
from the low-level, rapid variability of the He II 
A 4686 emission line in the spectrum of ^Puppis 
dEversberg et al.ll 1998D - Indirect evidence stems from 
the analysis of the spectra of O-type stars. Indeed, 
with observations from the FUSE satellite, it became 
clear that the modelling of lines of ions such as O V, 
N IV and P V with model atmosphere codes such as 
CMFGEN , is impossib le if clumping is not included 
(e.g. iBouret et all2005l) . 

The existence of such small-scale wind structures 
has important consequences on many diagnostics of 
mass-loss rates. Indeed, strong clumping often im- 
plies a substantial reduction, by factors between a few 
and about ten, of the overall mass-loss rate compared 
to a homogeneous stellar wind. In this context, our 
ignorance of the properties of these wind structures 
translates into uncomfortably large uncertainties on 
the actual mass-loss rates. For instance, whilst the 
space between clumps is generall y supposed to be to- 
tally empty, some recent studies (jZsargo et aL 2008; 
Sundqvist. Puis & Feldmeier 1 120101) showed that it is 
necessary to include the effects of a tenuous inter- 
clump medium to correctly predict the strength of var- 
ious lines in the observed UV spectra of O-type stars. 
By doing so, these authors even found that the ob- 
served profiles might again be compatible with the 
conventional mass-loss rates. 

Mass-loss rate diagnostics in different wavelength 
domains are strongly sensitive to different aspects of 
the wind fragmentation. For instance, H„ and radio 
emission are collisional processes which scale with 
density squared and are thus sensitive to the clumping 
factor regardless of the size of the clumps. In clumped 
winds with a void inter-clump medium, the velocity 
law is intermittent and some velocities are not present 
along a given sight line. T his porosity in the veloc- 
ity field is called 'vorosity' (Owocki 2011). The latter 
strongly modifies the Sobolev length and its impact is 
mainly visible on the UV resonance lines. 

The present paper deals with the impact of wind 
structure on the morphology of X-ray spectra of mas- 
sive stars. High-resolution Chandra and XMM - 
Newton X-ray spectra of presumably single O-type 
stars revealed less absorption of the X-rays by the 
cool wind than expected for a homogeneous wind 
with conventional mass-loss rates. This property can 
be explained in the case of fragmented stellar winds 



(e.g. iMiller et al] 120021: iKramer. Cohen & Owocki 



2003 



Oskinova, Feldmeier & Hamann 



2006t 



see 

also Giidel & Nazej 120091 for a review), although 
Cohe n et al.l (2010) argue in favour of an overall re- 
duction of the mass-loss rate of the wind without the 
need to include porosity effects. 

In a fragmented wind with optically thick clumps, 
X-ray lines are sensitive to porosity, which depends on 
clump continuum optical depth, and thus clump size 
as well as the clumping factor. That is why high- 
resolution X-ray spectra are important tools to study 
the structure of the wind and help us determine mass 
loss rates which must be consistent with the other 
wavelength domain analyses. Previous works in the 



X-ray domain (e.g. lOskinova. Feldmeier & Hamann 



2006), were performed on a line-by-line basis; the 
spectra were never analysed as a whole at once. Using 
multi-temperature plasma models in future work, the 
global fit of observations will give information con- 
cerning the different shock temperatures as well as the 
size and the position of the different X-ray emitting 
shells. 

With this in mind, we have designed a new code 
for the purpose of simulating the full high-resolution 
X-ray spectrum of a single massive star. Our method 
properly accounts for the effect of line blending on the 
line profiles, and allows testing the consistency of cur- 
rent models for the X-ray emission of massive stars 
over the full range of the X-ray domain. Therefore, 
our code is meant as a further step towards extracting 
the maximum amount of information from the high- 
resolution X-ray spectra. In this paper, we introduce 
the different free parameters and investigate their im- 
pact on single temperature synthetic X-ray spectra. We 
further discuss the main assumptions and physics in 
our code and compare different formulations of the 
wind fragmentation. In particular, we investigate for 
the first time the impact of a tenuous, but non-void, 
inter-clump medium on the synthetic X-ray spectra. 
Finally, we present the first application of our code to 
fit the region around the O VIII A 18. 97 A line in the 
XMM spectrum of ^Puppis and illustrate the degen- 
eracies on the stellar wind parameters that exist when 
a single line is fitted. In future work, we will apply this 
code to the full X-ray band spectra of ^Puppis. 

2. Assumptions and main concepts of the code 

Before using a multi temperature plasma model to 
analyse the X-ray spectra of O-type stars in detail, we 
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need to calculate single temperature synthetic spec- 
tra with the most recent stellar wind theories (emis- 
sion and absorption) and with as few as possible free 
parameters. We decide to follow the wind embed- 
ded shocks scenario. Consequently, we need to know 
where the shocks arise in the wind and what is the tem- 
perature of the resulting post-shock hot plasma. Us- 
ing up-to-date emissivities, we can then determine the 
X-ray emission of a hot plasma as a function of tem- 
perature, abundances and quantity of matter present in 
the shocks. Next, we determine how the cool wind 
material absorbs the X-ray photons and how the wind 
structure (i.e homogeneous or porous) impacts on the 
emergent X-ray flux. In the following subsections, we 
describe the tools and assumptions considered to build 
our code and to calculate synthetic spectra. 

2.1. The plasma emission spectrum 

Previous works on X-ray line profile fitting have al- 
ways normalized the line profiles, i.e the overall line 
flux was never physically modelled. We attempt a dif- 
ferent approach here, by using the true X-ray emissiv- 
ity of an optically thin plasma. For this purpose, we 
use AtomDlf] which models the emission from a col- 
lisionally ionized, optically thin isothermal plasma in 
thermal equilibrium. Many processes are included in 
the model: collisional excitation, ionization, and di- 
electronic recombination drive the ionization balance 
and emission spectra from the plasma. Continuum pro- 
cesses such as two-photon decay and bremsstrahlung 
are also included. Photon-driven processes such as 
photo-excitation and photo-ionization are not included 
in the model. The assumption of an optically thin 
plasma requires a low density (< 10 13 cirT 3 ), there- 
fore three body effects are not included. The standard 
AtomDB model data has been recalculated on a finer 
0.025 keV/step electron temperature grid for this anal- 
ysis and we use an energy resolution of 0.2 eV. 

In addition to the strongest lines that dominate the 
spectrum, the high-resolution theoretical spectra re- 
veal the presence of numerous weak lines and also 
discontinuities in the continuum (Fig. [TJ. Their con- 
tribution to the emergent flux is non-null and it is nec- 
essary to take them into account for the computation 
of emergent spectra as they certainly modify the emis- 



2 This was previously known as the APEC model ISmith et alj 
l200ll) . The interested reader is invited to visit the section 'physics' 
on the AtomDB website for a precise description of the calcu- 
lation of the continuum and lines in a hot plasma at the url 
http://www.atomdb.org/ 



sivities and profiles of the strongest observed lines. 
When fitting actual observations, one further needs 
to account for the presence of plasma at differ- 
ent temperatures, located at different places in the 
wind. Therefore the normalization of the model 
will result from the combination of several model 
spectra weighted according to their relative emis- 
sion measure. This will be dealt with in our forth- 
coming paper on model fits of ^Puppis spectrum. 
In addition, in the X-ray emitting plasma, some ions 
are present in a helium-like form. The treatment 
of the transition between the 2 3 S 1 and 2 3 P, levels 
must be done separately. Indeed, the presence of a 
strong UV radiation modifies the populations of the 
metastable upper level of the forbidden line via a 
pum ping to the upper level of the inter-combin ation 
line felument hal et alJll972t IPorquet et alj|200ll) . We 



will discuss this fundamental effect in our more com- 
plete analysis of ^Puppis presented in a forthcoming 
paper (Herve et al., in preparation). 

We consider the hot plasma to be embedded in the 
cool wind and to move outwards along with the latter. 
We decompose the emitting volume into small cells 
where we calculate the X-ray emission of each cell and 
its shift in wavelength in the observer's frame of refer- 
ence due to the overall wind motion. 

2.2. Wind optical depth and mass-absorption co- 
efficient 

The X-ray emitting plasma is assumed to be op- 
tically thin to X-rays. An X-ray photon emitted by 
the hot plasma, either escapes from the wind or is ab- 
sorbed by the cool wind material. Thus the radiative 
transfer is reduced to the calculation of t, the opti- 
cal depth along the line-of-sight du e to the cool wind 
materi al. Following the approach of lOwocki & Cohenl 
the absorption by the stellar wind is essen- 
tially characterized by the t* parameter and the mass- 
absorption coefficient k, which are linked by the rela- 



tion t» = 



kM 



These absorption coefficients are de- 
pendent on the wavelength and the abundances of the 
different elements. To evaluate the mass-absorption 
coefficients, we use the radiative transfer code CM- 
FGEN (Hillier & Miller 1998) which computes the 
ionization structure of the cool wind as a function of 
temperature, abundances of many elements and mass- 
loss rate. In the ionization calculations, CMFGEN 
accounts for different physical processes such as col- 
lisional excitation, photo-ionisation by photospheric 
light and X-rays, and radiative and dielectronic recom- 
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bination. From the known ionization structure, we can 
then compute k. In the X-ray domain, the opacity is 
dominated by the photo-ionization of the metals. This 
absorption is by far the dominant source of opacity and 
we can safely neglect other physical processes such as 
free electron scattering. 

The wavelength dependence of k (Fig. |2]i and there- 
fore r», is quite sensitive to the model elemental abun- 
daces and to the ionization fraction of singly ionized 
helium. In this respect, we stress that the large values 
of the mass absorption coefficient, used in the present 
paper, reflect the fact that the sum of the CNO abun- 
dances in the best-fit model of the UV and optical spec- 
trum of ^Puppis is supersolar (Bouret et al., in prepa- 
ration). 

The detailed behaviour of a: as a function of radius 
and wavelength (especially at longer wavelengths) also 
depends upon the ionization structure of the wind and 
more specifically on the details of the rec ombination 
of He 2+ into He + ( Leutenegger et al]|2010b . However, 
for simplicity, we assume here that k is independent of 
the position in the winc0 and we evaluate its value at 
r = 10R*. In our detailed analysis of the RGS spec- 
trum of ^Puppis (Herve et al., in preparation), we will 
drop this assumption and rather introduce a radial de- 
pendence of K. 

Finally, the values of the mass-absorption coeffi- 
cient or the wind optical depth as a function of the 
wavelength are used in Eq.|3]and|4](see Sect. |2.3| be- 
low) to determine the absorption of X-rays by the stel- 
lar wind. 

2.3. X-ray absorption by a fragmented wind 

Previous studies of stellar winds revealed that 
they are not homogeneous but fragmented. The ab- 
sorption of X-rays by a clumpy wind has been ad- 
dressed by several authors who adopted slightly dif- 
ferent formalisms for the same phenomenon. The 
porosity length prescription is based on the con- 
cept of the mean free path of photons between two 
successive interacti ons with clumps, and the opac- 
ity of the latter dOwocki & Cohenl Eooit) . In the 
fragmentation frequency prescription, the opacity of 
the wind in the X-ray band is estimated by the fre- 
quency of clumps passing through a ref erence radius 
(Osk inova. Feldmei er & Ha mann I l2006l) . In this sec- 
tion, we briefly summarize the main aspects of these 



models. 

At the microscopic level, the absorption of X-rays 
by a stellar wind is mainly due to the bound-free tran- 
sitions from the K and L energy levels. To start, 
we assume that the clumps are spherical and very 
small. We consider clumps over a wide range of 
optical depth, from optically thin to optically thick. 
In a situation where the individual clumps are opti- 
cally thick, the effective opacity is essentially deter- 
mined by the geometrical properties of the clumps 
( Feld meier. Oskinova. & Ham ann 120031) and can be 
written K e ff = — = — where k is the mass-absorption 
coefficient, I is the geometrical scale of the clump, m c 
its mass a nd r r its optical depth. This latter can be 
written as (Owocki & Cohenll2006l) : 



/ kM \ I 
Tc=K(p) f = 4^v7 



(1) 



where (p) is the mean density of the stellar wind, / 
the volume filling factor, M is the mass-loss rate of the 
star and v the velocity of the wind at the radius r which 
is determined by a /Haw (i.e v(r) = v^il - —)P with 
ft = 1 throughout this work). 

Under this assumption, it appears that the effective 
reduction depends on the ratio between the clump 
scale and filling factor, a lso called the porosity length 
(lOwocki & Cohen 2006). The concept of effective 
opacity can be extended t o the optically thin limit us- 
ing a simple bridging law ( Owocki & Cohen| [2006): 



K eff 
K 



1 + Tr 



(2) 



Using this result in a steady state (i.e. non variable) 
wind, one obtains the wind optical depth in co nven- 
tional (p, z) coordinates (lOwocki & Cohenl 2006): 



r(p,z) 



7<r~ 



R»dz' 



R t ) + T t h{r') 



(3) 



where r' — p 2 + z' 2 



T» = 



kM 



is the optical depth 



This assumption allows us to take k out of the optical depth integral 
(Eq.fJJ and thus to achieve an analytical evaluation of this integral. 



parameter of the wind, Vaa is the terminal velocity of 
the wind, R* is the radius of the star and h(r) = 4 the 
porosity length. In the following we adopt h(r) — hxr. 

The fragmentation of stellar wind shells is a conse- 
quence of hydrodynamical instabilities. Rather than 
being spherically symmetric, the resulting clumps 
could actually have a flattened, elongated shape. Using 
the number «o of fragments passing through some ref- 
erence radius per unit time, also called the fragmenta- 
tion frequency, and the opacity of each fragment at this 
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referen ce radius, lOskinova. Feldmeier & Hamann 
(|2006l) showed that the wind opacity can be written 
in the general case as: 



, z) = «o J 



{\-e- T ')W)\ 



dz' 
v(r') 



(4) 



where r ; is the optical depth of a flattened fragment j 



along the line-of-sight and p(r') 



represents 



its orientation. The optical depth along the line-of- 
sight: 



Ti = 



\p\ 4tt r 2 n 



(5) 



with k the mass-absorption coefficient and t . the av- 
erage r adial optical depth of a fragment located at dis- 
tance r (lOskinova. Feldmeier & Hamannl2004l) . In the 
case of an isotropic opacity (i.e. small-scale spheri- 
cally symmetric clumps), p is taken equal to unity (i.e. 
the line-of-sight crosses each clump along a radial di- 
rection). 

In Section 3.1, we compare the profiles obtained 
with t given by Eq.|4] to those computed with Eq.[3] 
first in the case of an isotropic opacity where p is taken 
equal to one, then in a case of an anisotropic opac- 
ity where the clumps are assumed to be flattened (i.e 

2.3.1. Differences and common features of the differ- 
ent prescriptions 

As we will show below, the porosity length pre- 
scription and the fragmentation frequency prescription 
are essentially equi valent to first order, al though in 
their original work, lOwocki & Cohen (2006) concen- 
trated on spherical clumps, whilst lOskinova. Feldmeier & 
(2006) accounted for flattened structures by introduc- 
ing the line-of-sight orientation via the u parame tefl 
Actua lly, the geometry assumed by lOwocki & Cohenl 
(2006) is even more idealized than a pure spherical 
geometry since it assumes that the line-of-sight al- 
ways crosses the clump along the radial direction. 
Strictly speaking, this is valid only if the clumps are 
very small. The main differences between the two 
approaches hence concern the geometry of the wind 
fragments and the way the optical depth of a clump 
is included in the model, i.e. Eq.|2] for the porosity 
prescription, and the (1 - e~ T ') term in Eq.|4]for the 
fragmentation frequency. For small optical depths and 



isotropic clumps (jj, = 1), both models should be es- 
sentially equivalent. 

In a first step, we reproduce the profile of various 
lines (Table[T|i for each prescription and compare them 
with each other and with the results for a homogeneous 
wind model. We consider a distribution of spherically 
symmetric clumps and we thus set p = 1 in Eq.0 and 
[5] Our general results agree well with those of previ- 
ous studies. First, the larger the porosity of the wind, 
the larger the mean free path of a p hoton in the stellar 
wind and the lower the absorption dOwocki & Cohen 



2006). Second, the larger the fragmentation frequency, 
the lower the opacity in the fragmenta t ion pr escription, 



see 



Oskinova, Feldmeier & Hamann (2006). 



We have simulated four specific lines with np = 1.7 
10~ 4 s " 1 adopted from lOskinova. Feldmeier & Hamann 
d2006h . This parameter is found by the authors to 
best fit the O VIII A 18.97 A line in the HETG and 
RGS spectra of ^Puppis. Then we derive the porosity 
length, h, needed to reproduce the line profiles com- 
puted with this fragmentation frequency (Fig. |4j. For 
different values of t», hence different values of wave- 
length (see Table 1), the two line profiles are, in each 
case, very close. Note however that we obtain dif- 
ferent values for the porosity length when matching 
different lines computed with the same fragmentation 
frequency: we have to decrease the porosity parameter 
when t* increases. In other words, when we com- 
pare simulated global spectra using a single porosity 
parameter to those obtained for a specific value of 
the fragmentation frequency, we find that the spectra 
are exactly the same in the wavelength band near the 
line used to find the parameter h that best matches 
the value of «o (in our case the O VIII A 18.9 A 



4 Note that this parameter could also easily be inserted into Eq.f2] 



Hamatane), but differ over the other wavelength domains 
(Fig-Ell. More specifically, the porosity model predicts 
a slightly lower flux in the short wavelength domain 
and a higher flux in the long wavelength band. 

There are some underlying assumptions in both 
models, and with current data it seems difficult to tell 
which prescription provides the better description of 
reality. The porosity prescription has the advantage to 
provide an analytical form of the opacity that speeds 
up the computation. Therefore, in the following, we 
concentrate mainly on the porosity prescription. The 
overall results and conclusions obtained in the study 
for the different parameters, in the next section, are 
the same for the fragmentation model in an isotropic 
(p = 1) case. 
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3. The impact of model parameters on the emer- 
gent spectrum 

For each simulation in this section, we use a termi- 
nal velocity Voo of 2250km s _1 . We further assume a 
kl plasma of 0.217keV, solar abundances for the chem- 
ical composition of the plasma, a stellar radius R* 
of 18.6 7? and a mass-loss rate of 4.2 10 _6 M o yr _1 . 
These stellar and wind parameters correspond to those 
used by Oskinova et al (2006) in their work to fit lines 
of £ Puppis, except for the plasma temperature which 
corres ponds to the closes t value in AtomDB to the re- 
sult of Iffillier et ail d 19931) . 

In the remainder of this section, we simulate a series 
of spectra for such a single O star in order to illustrate 
the sensitivity of the synthetic spectra on the different 
model parameters. 

3.1. Porosity 

The first explored parameter is the porosity length. 
The impact of this parameter for broad band spectra 
was never presented, in the literature. Fig.|6]illustrates 
our results. We have chosen to explore a wide range 
of values of the parameter h. As expected, we find 
a strong impact on the emergent flux. At the short- 
est wavelengths, where k is lowest, the porosity length 
has little impact on the spectrum. However, the situa- 
tion is radically different at longer wavelengths, where 
the wind opacity k becomes huge (Figf3] Tab[TJ. Here, 
the impact of the porosity length is much more impor- 
tant. We find differences of several orders of magni- 
tude between the spectrum for a homogeneous wind 
and the simulations with the largest porosity. For a 
highly porous wind, the variations of h have less im- 
pact as the wind is already so fragmented that the pho- 
tons have few remaining interactions with the clumps. 

As a second step, we consider the effect of the shape 
of fragmented shells. For this purpose, the factor // 
in Eq.[4]and[5] which represents the orientation of the 
fragmented shells with respect to the line-of-sight, is 
re-introduced into our simulations. Compared to the 
isotropic case, the opacity of the wind decreases (see 
Fig-0- The radiation from the back side of the star 
is less attenuated than in the case of isotropic clumps. 
Indeed, the fact that fragments are flattened and that 
photons coming from the rear side of the wind cross 
these shells mainly perpendicularly implies that they 
have to cross less matter, and are thus less absorbed, 
than in the isotropic case. As a result, the total flux 
and intensity are higher in the anisotropic case. Also, 



the profiles are less asymmetrical and the full width at 
half maximum increases. The consequences are more 
intense spectra (by a factor 1 .2 in the short wavelengths 
and 3.0 in the longer wavelengths domain) and an in- 
creased overlap of lines in the anisotropic model (see 
Fig.®. 

3.2. Outer radius of the emitting region 

Another parameter of our model is the size of the 
X-ray emitting plasma (R out )- As the emissivity is pro- 
portional to p 2 , the outer parts of the wind contribute 
less to the total intrinsic (i.e. before wind absorption 
is accounted for) X-ray emission, but cannot be ne- 
glected for the wind-absorbed emission as can be seen 
for specific line profiles in Fig. [9] The effect is by far 
strongest in the long wavelength domain where /q is 
largest. The intensity of the lines in Fig. [9] varies by 
a factor between 1 . 1 and ~ 80 (from shorter to longer 
wavelengths) when the outer radius of the emitting vol- 
ume is multiplied by a factor 2.0. Moreover, extend- 
ing the emitting region to larger distances implies that 
the velocity of the X-ray emitting plasma (assumed to 
move along with the cool wind) becomes larger if the 
outer radius is still located inside the wind acceleration 
zone. This in turn increases the full width at half max- 
imum of the lines as well as the line width at the con- 
tinuum level. The effect is best seen in the red wing of 
the line which is most strongly depressed by the wind 
absorption. With a larger emitting volume, a larger 
fraction of the redshifted photons from the rear outer 
parts of the wind escape the wind since they cross lay- 
ers of lower density. The impact of the increase of the 
size of the emitting region on the overall spectrum can 
be seen in Fig.fTOl 

3.3. Inner emission boundary 

The location of the inner boundary of the emission 
region is one of the open questions in the understand- 
ing of the X-ray emission of massive stars. However, 
a variation of this parameter can have important con- 
sequences. As we pointed out in the previous section, 
the density of matter increases as one moves inwards 
in the stellar wind, and therefore the contribution to 
the intrinsic X-ray emission increases as well. Still, 
as far as the observable emission is concerned, part 
of this increase in emission measure is compensated 
by the increase of the amount of absorbing material 
that the photons have to cross before they escape. If 
we adopt an inner radius of 2.5/?* instead of 1.5/?*, 
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with an outer radius of 10. for the emitting plasma, 
we note that the line intensity in the short wavelength 
range, where t» is low, decreases roughly by a factor 
two. In the long wavelength domain, where r» is high, 
the reduction of the stellar wind absorption is more 
important than that of the emission. In consequence, 
the inner part of an X-ray emission shell is fully ab- 
sorbed by the cool material. The result is that the lines 
at longer wavelength are rather insensitive to the lo- 
cation of the inner boundary of the emission region 
(FigfTTT). The same qualitative considerations hold for 
the overall spectrum which is shown in Fig.fT2l 

We note that the FIR triplet of helium-like ions 
provides in principle a po werful diagnostic of th e 



inner emission boundary dLeutenegger et al.l |2006). 



However, in the case where the observed X-ray emis- 
sion stems from a multi-temperature plasma with each 
component possibly having its own value of the inner 
boundary, it must be stressed that the resultant triplets 
are the sum (weighted by the emission measures) of 
the triplets in the various components. This situa- 
tion renders the interpretation of the FIR diagnostics 
more complex than expected from a single temperature 
plasma. 

3.4. The contribution of the inter-clump medium 



In their work, IZsargo et al.l d2008) have shown the 
need to introduce a homogeneous inter-clump medium 
to better fit the profile of the O VI /LI 1032,1038 
doublet with the model atmosphere code CMFGEN. 
These authors conclude that the inclusion of a small 
percentage (~ 5%) of homogeneous material in be- 
tween the clumps is needed to achieve a good fit 
of the line. Similar conclusions were reached by 
Sundqvist. Puis & FeldmeTerl d2010h . 

What is the consequence of such a tenuous, but 
homogeneous, inter-clump medium on the X-ray line 
profiles? To answer this question, we incorporate this 
situation in our model. In practice, we work with 
a constant mass-loss rate and we consider that the 
size of the clumps is small. We distribute an im- 
portant percentage of matter (90%, 95% or 97%) in 
clumps and we compute the optical depth due to a to- 
tally clumped wind. As a first assumption, we con- 
sider that the ionization is the same in the cool clump 
as in the inter-clump medium^ Next, we compute 



the optical depth for an homogeneous wind with the 
remaining matter. Finally we take the sum of the 
two contributions and determine the absorption for a 
mixed wind model. As a homogeneous wind absorbs 
more photons, the intensity of the line profiles of the 
mixed model should actually be lower than for a fully 
clumped model (see Fig.[T3l and the impact is stronger 
at longer wavelengths, and hence at larger opacities 
k. These considerations are confirmed by our sim- 
ulations of the overall X-ray spectrum (see Fig.fPfl). 
As a second step, we attempt to fit the spectrum of the 
mixed model with 97% of matter clumped to a model 
spectrum for a fully clumped wind. Fig. [TSI illustrates 
the results: the individual line profiles of the mixed 
model are almost perfectly fitted by a fully clumped 
model with a slightly lower porosity parameter. This 
has important consequences in practice: whilst the 
X-ray line profiles might allow us to infer the pres- 
ence of clumps and to constrain their geometry, it 
seems almost impossible to distinguish between a fully 
clumped wind and a mixed model without constraints 
from observations at other wavelengths (e.g. UV, op- 
tical). We thus conclude that X-ray line profiles alone 
cannot reveal the presence of an inter-clump medium. 

4. First application to the O VIII A 18.97 A line of 
£Puppis 

In the previous sections, we have discussed the im- 
pact of different parameters and of different fragmen- 
tation prescriptions on the line profiles and on the 
emergent X-ray spectra. In this section, we use our 
code to fit the O VIII Lyman a doublet ,118.97 A and 
the neighbouring N VII A 19.36 A, A 19.82 A lines of 
the massive star £Puppis and determine some stellar 
wind parameters. We decided to work on an extended 
region around the strong O VIII line since Fig. [T] re- 
veals the presence of weak lines that affect the wings 
of the O VIII line and produce a pseudo-continuum 
that contributes of order 5% of the flux of the O VIII 
line. For this purpose, we have extracted, reduced 
and combined 18 RGS spectra of fPuppis from the 
XMM - Newton archive (Naze . Flores & Rauw hoi a 
the different aspects of the data reduction are treated 
in detail in that paper). We also apply the absorp 



,-2 



5 Constraining the ionization of the inter-clump medium requires hy- 
drodynamical simulations that are beyond the scop e of the present 
paper (see also Sundqvist, Puis & Feldmeier 2010 and Zsargo et al. 



tion by the interstellar medium (Nh = 8.9 x 10 cm 
iDiplas & Savage 1 19941) to the synthetic spectra. Then, 
we convolve the spectra with a Lorentzian profile with 
a half-width at half-maximum of 0.07 A to account for 
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the instrumental resolution of the RGS spectrograph 
onboard XMM - Newton. Finally, we compare the re- 
sults to the observation. 

As shown in Fig. [16] and [17] we obtain a degeneracy 
of the stellar wind parameters. Indeed, for the O VIII 
line profile, we demonstrate that we can use a homo- 
geneous or a porous wind model to fit the data equally 
well (see Fig. [16] and Table 2) though with different 
values for the hot gas filling factors. Note that these 



values are in agreement with the results o flHillier et al. 
dl993h . We only need to adapt the mass-loss rate and 
the abundance of oxygen to obtain similar results. In 
a denser wind (i.e. with a larger mass loss rate), the 
increase of absorption dominates the increase of X- 
ray photon production at longer wavelengths where the 
opacity is largest. In consequence we have two ways 
to increase the number of X-ray photons which escape 
from a denser wind. The first solution is to increase 
the number of X-rays resulting from shocks. To first 
order, we could imagine that the number of shocks, 
and hence the hot gas filling factor, would be larger 
in a denser wind (Tab.[2] models A and C). Alterna- 
tively, we can increase the number of photon which 
escape. Consequently, we have to increase the free 
mean path of photons in the stellar wind. This can be 
achieved with a more porous wind (Tab.[2] models A 
and D). Normally, the line profiles of the two previous 
solutions should be different. Unfortunately, the in- 
strument resolution is not sufficient to distinguish the 
subtle differences. 

Furthermore, we can also fit this line equally well 
with a model where the inter-clump medium is not 
void (Fig.Ha on the left). 

It is thus obvious that a single line does not per- 
mit to discriminate between models. These prelimi- 
nary results therefore emphasize the necessity to work 
on the full X-ray spectrum. The analysis of the full 
RGS spectrum is beyond the scope of the present pa- 
per, but will be presented in a future work (Herve et al., 
in preparation). As we demonstrated in Sect. 3, each 
parameter has a different impact over different parts of 
the spectrum and we are therefore confident that some 
degeneracies on e.g. porosity or mass-loss rate, will 
then be lifted. 

Nevertheless, already on this small region of the 
spectrum, three results are found. First, the possibil- 
ity of an anisotropy (i.e., fi + 1, in the fragmenta- 
tion prescription) of the wind absorption appears un- 
likely. Indeed, the flux from the back side of the 
star would be too large (Fig. [7] on the top right) and 



hence the width of the synthetic line emission would 
be too large in comparison to the observations (Fig. [18] 
on the right). Second, a good fit to the O VIII line 
and its close neighbours is only found if one uses an 
overabundance of N and a depletion of O (Tab. [2] and 
Fig.[T6l). Third, a quick look at a wider wavelength 
band (FigJT9b shows that a single temperature plasma 
cannot reproduce the whole spectrum of ^Puppis. A 
multi-temperature plasma model is required for a full 
analysis of the X-ray spectrum. Thus the observed line 
profiles and intensities are a combination of the contri- 
bution of the different temperature plasmas present in 
the wind. 

Finally, it should be noted that part of the degen- 
eracy of the results is due to the limited resolution of 
present-day high-resolution X-ray spectrographs. In 
Sect. 3, we show the impact on a line profile of a 
change in porosity length or inner radius for example, 
but the convolution by the response matrix of current 
instruments washes out part of the difference induced 
by the different parameters. 

5. Conclusions 

In this paper, we have studied the impact of the 
two different prescriptions of stellar wind fragmenta- 
tion on the emergent broad-band X-ray spectrum. For 
spherical clumps with an isotropic opacity, the dif- 
ferences between the two prescriptions are the most 
obvious when comparing the overall X-ray spectrum. 
Indeed, in this case, we find that a given simulation 
with a single value of the fragmentation frequency 
cannot be represented by a spectrum calculated with 
a single value of the porosity parameter. However, in 
practice, with the quality of current high-resolution 
X-ray spectra of O-type stars, it seems almost im- 
possible to distinguish the two models in a real case 
when working on a line-by-line basis. This is because 
the differences are often small and are washed out 
by the coarse resolution of present-day instruments. 
Moreover, the differences are most easily seen when 
comparing the shape of the full spectrum which is sel- 
dom done. In addition, one must keep in mind that 
the shape of spectrum is also affected by the amount 
of interstellar absorption and the fact that the X-ray 
emission from a real stellar wind most likely spans 
a range of plasma temperatures rath er than a single 
temperature as assu med here (see e.g. lZhekov & Palla 



2007; iNaze 1 120091) . In both prescriptions we find 



that the X-ray emission is more intense in a frag- 
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mented stellar wind than in a homogeneous wind. 
The introduction of an anisotropic opacity in the 
fragmentation frequency prescription yields spectra 
that are more intense than in the isotropic case and 
of course than in an homogeneous wind. Sophis- 
ticated, 3-dimensional hydrodynamical calculations 
might eventually be needed to clarify the issue of the 
geometry of the clumps. 

As pointed out previously, both models rely on 
a series of assumptions and, currently, we do not 
have enough knowledge on the stellar wind frag- 
mentation to privilege one model over the other. 
Our investigation of the impact of the various model 
parameters on the emerging spectrum showed that both 
the position and the size of the X-ray emitting plasma 
are very important. An increase in the size leads to 
a more luminous spectrum, although we have shown 
that the amplitude of the intensity increase depends 
upon wavelength. 

For a given mass-loss rate, the differences be- 
tween fully clumped and homogeneous wind mod- 
els are quite large. However, in light of the sim- 
ulations presented here, it will be very difficult 
to estimate correctly the value of the porosity if 
we include the effect of an inter-clump medium. 
In fact, the incorporation of a few percent of the 
wind material in the inter-clump medium yields the 
same emergent spectrum as obtained for a purely 
clumped model with a slightly less fragmented wind. 

Finally, we present preliminary fits of the region 
near the O VIII A 18. 97 A emission line of ^Puppis 
with our code. Our first results reveal that anisotropic 
clumps are unlikely, and that non-solar abundances 
(overabundance of N, depletion of O) are needed to 
reproduce the X-ray spectrum. However, there is also 
a degeneracy of stellar wind parameters. Even though 
a porous wind is closest to the observational data, we 
can not totally exclude an homogeneous wind. The 
mass-loss rate, the porosity parameter and the hot gas 
filling factor are strongly tied. However, simultaneous 
fits of the total spectrum will likely allow to reduce the 
degeneracy on the porosity and other parameters. But, 
it appears that a single temperature X-ray plasma is in- 
sufficient to fit the global spectrum of ^Puppis. This 
result implies the necessity to use several temperatures 
and we will have to include the contribution of each 
temperature plasma to the line intensity and profile of 
O VIII A 18. 97 A (and all the other lines) but also to the 
continuum. 



Table 1 : Lines considered in this study and the wind 
optical depth extracted from Fig. [3] and corresponding 
to our stellar and wind parameters 



Line 


Wavelength (A) 


T* 


Si XIV 


6.18 


1.8 


O VIII 


18.97 


21.9 


C VI 


33.73 


46.2 


Si VIII 


61.02 


187 



1 We are aware that the Si VIII A 61.02 line will 
be impossible to observe in practice because it 
is expected to be intrinsically faint and most of 
all because of the huge absorption by the inter- 
stellar medium at these wavelengths. However, 
here we use this line to illustrate the behaviour 
of lines at high r, values. 




Wavelength in Angstrc 



Fig. 1 . — Theoretical emission of a plasma with a tem- 
perature of kT=0.273 keV and with a ^ = 0.3 and 

jf- — 4. by number. 
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6.14 6.16 6.18 6.20 6.22 
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18.85 18.90 18.95 19.00 19.05 19.10 
Wavelength in Angstrom 



33.5 33.6 33.7 33.8 33.9 
Wavelength in Angstrom 




}.6 60.7 60.8 60.9 61.0 
Wavelength in Angstrom 



Fig. 2. — Evolution of the mass-absorption coefficient 
as a function of the wavelength. In black, the mass 
absorption coefficient with our preliminary results on 
the abundances determined in the X-ray spectrum of 
^Puppis reveals a less opaque wind than with the 
abundances determined in the most recent UV/optical 
analysis (in red, J-C Bouret, private communication, 
Bouret et al., in preparation). 



Fig. 4. — Illustration of the impact of porosity (h = 
0.07, 0.07, 0.055 and 0.043 R„ dashed red line, from 
left to right and top to bottom) on the profiles of four 
specific lines (see Table 1). The profiles are compared 
with those computed in the fragmentation prescription 
for no = 1.710~ 4 s (black solid line) and for a homo- 
geneous wind model (dot-dashed green line). 



250 




Fig. 3. — r* computed with solar abundances, M = 
4.2 lO^Moyr 1 , v M = 2250^8-' andR, = 18.6 R . 
These values are used in Section 3. 




10 15 20 25 30 35 40 

Wavelength in Angstrom 



Fig. 5. — Comparison between a synthetic spec- 
trum computed with the porosity prescription adopt- 
ing h = 0.07 R* (dashed red line), a spectrum com- 
puted with the fragmentation frequency theory with 
no = 1.710~ 4 s _1 (black solid line) and a homogeneous 
wind (dot-dashed green line). 
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Fig. 6. — Effect of porosity on the emergent spectrum. 
In long-dashed magenta line, the results for a homo- 
geneous wind, in dashed-double-dotblue line, dashed- 
dot green line, dashed red line and solid black line the 
spectra for different porosity parameters (0.02, 0.07, 
0.1 and 0.15/?*, respectively). A porosity of 0.15/?* 
corresponds to the brightest model, whilst 0.02/?* cor- 
responds to the most absorbed (hence weakest) model. 
See the online version of the article for a colour figure. 



o.Qoi L i i i i i i i i i i i i i i i i i i i i i i i i i i i i i l>, j\JV| V 1 

10 15 20 25 30 35 40 

Wavelength in Angstrom 

Fig. 8. — Same as Fig. [7] but this time for the entire 
X-ray spectrum. 




33.6 33.7 33. J 
Wavelength in Angstrom 



Fig. 7. — Impact of the anisotropy of the opacity (frag- 
mented shells seen under a direction cosine ft) on the 
emergent line profile, shown as the dashed red line, for 
specific lines (see Table 1). The results are compared 
to the line profiles computed for an isotropic opacity 
(/i = 1; solid black line) withno = 1.710 _4 s _1 . 



Fig. 9. — Impact of the size of the emitting plasma 
on the line profiles at different wavelengths. The solid 
black (respectively dashed red) profiles correspond to 
an outer radius of the emitting region of 10/?* (resp. 
5/?*). 
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Fig. 10. — Impact of the size of the emitting plasma on 
the overall X-ray spectrum. The colours have the same 
meaning as in Fig. [9] 
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Fig. 1 1 . — Importance of the inner boundary of the 
X-ray emission region. The dashed red (respectively 
solid black) line profiles illustrate the results of our 
simulations for different lines assuming ro = 2.5 R* 
(resp. 1.5/?,) and on outer radius of the emitting re- 
gion of 10/?*. The huge t, value at large wavelengths 
(Fig|3] Tab[TJ makes the wind opaque to X-ray pho- 
tons. Consequently the X-ray photons produced in the 
inner part of the emitting shell are all absorbed by the 
cool wind material (bottom panels and right top panel). 
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Fig. 12. — Same as Fig.[TT] but this time for the overall 
X-ray spectrum. 
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Fig. 14. — Same as Fig. [13] but illustrating the overall 
spectrum between 20 and 35 A. 
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Fig. 13. — Impact of the presence of an inter-clump 
medium. The solid black line profiles correspond to 
a fully porous wind (with h - 0.10/?*), whilst the 
dashed red profiles yield the results for a fully homoge- 
neous wind. The dot dot dashed blue, dot dashed green 
and long dashed magenta profiles display the situation 
for 3 different mixed models with respectively 90%, 
95% and 97% of the wind material in clumps. 
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Fig. 15. — Comparison between a mixed model with a 
fully clumped model (see text for details). 



13 



This research is supported by the FNRS (Bel- 
gium) and by the Communaute Franchise de Belgique 
- Action de recherche concertee (ARC) - Academie 
Wallonie-Europe and an XMM+INTEGRAL PRODEX 
contract. (Belspo) 

REFERENCES 

Blumenthal, G. R., Drake, G. W. F, & Tucker, W. H. 
1972, ApJ, 172, 205 

Bouret, J.-C, Lanz, T, & Hillier, D.J. 2005, A&A 438, 
301 

Cohen, D.H., Leutenegger, M.A., Wollman, E.E., 
Zsargo, J., Hillier, D.J., Townsend, R.H.D., & 
Owocki, S.P. 2010, MNRAS, 405, 2391 

Diplas, A., & Savage, B.D. 1994, ApJS, 93, 21 1 

Dopita, MA. 2008, in Massive Stars as Cosmic En- 
gines, Proc. IAU Symp., 250, 367 

Eversberg, T, Lepine, S., & Moffat, A.F.J. 1998, ApJ, 
494, 799 

Feldmeier, A., Oskinova, L. & Hamann, W.-R. 2003, 
A&A, 403,217 

Giidel, M., & Naze, Y. 2009, Astron. Astrophys. Rev., 
17, 309 

Hillier, D.J., Kudritzki, R.P, Pauldrach, A.W., Baade, 

D. , Cassinelli, J.P, Puis, J., & Schmitt, J.H.M.M. 
1993, A&A, 276, 117 

Hillier, D.J., & Miller, D.L. 1998, ApJ, 496, 407 

Hirschi, R. 2008, in Clumping in hot-star winds, eds. 
W.R. Hamann, A. Feldmeier, & L. Oskinova, 9, 
arXiv0709.0392 

Kramer, R.H., Cohen, D.H., & Owocki, S.P. 2003, 
ApJ, 592, 532 

Leutenegger, M.A., Paerels, F.B.S., Kahn, S.M., & Co- 
hen, D.H. 2006, ApJ, 650, 1096 

Leutenegger, M.A., Cohen, D.H., Zsargo, J., Martell, 

E. M., MacArthur, J.P, Owocki, S.P, Gagne, M., & 
Hillier, D.J. 2010, ApJ, 719, 1767 

Miller, N.A., Cassinelli, J.P, Waldron, W.L., MacFar- 
lane, J.J., & Cohen, D.H. 2002, ApJ, 577, 951 

Naze, Y. 2009, A&A 506, 1055 




0.0020 






I 0.0015 






"1 0.0010 






0.0005 







5 19.0 IS 

Wavelength in Angstrom 




Fig. 16. — Fit of the spectral region around the O VIII 
/U8.97 line of £Puppis. There exists a degeneracy on 
the determination of the mass-loss rate, the porosity 
of the wind and the hot gas filling factor. Top panels: 
Profiles derived for a homogeneous wind (on the left, 
model A in dashed red line) and a porous wind with the 
same mass-loss rate (on the right, model B in dashed 
red line) compared to the observational data (in solid 
black line). Bottom panels: same but for a different 
mass-loss rate (model C on the left, model D on the 
right). Model parameters are given in Table 2. 




Fig. 17. — Same as Fig. [16] but zooming in on the 
O VIII A 18. 97 line. 
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Table 2: Stellar and wind parameters for the fit of the O VIII /I18.97 A line in the RGS spectrum of ^Puppis with 
kT plasma = 0.217 keV, R t = 18.6 R Q , R inner = 1.5R t , R outer = 10/?*, v m = 2250kms _1 , £ = 0.3 (by number) and 

= 4.0 (by number). 



Model 


M 

(lO^Moyr 1 ) 


fhot gas 


h 
(*.) 


note 


2 

Xv 


A 


3.0 


0.033 


0.00 




4.10 


B 


3.0 


0.028 


0.10 




5.44 


C 


4.0 


0.036 


0.00 




4.39 


D 


4.0 


0.029 


0.10 




3.92 


E 


4.0 


0.028 


0.10 


5% of inter-clump 
matter 


5.31 


F 


4.0 


0.026 


n = 1.710-V 1 


H + 1 (fragmentation 
prescription) 


8.32 



0.0025 
0.0020 
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Fig. 18.— Left: Comparison of the O VIII ,118.97 A 
line profile of the non-void inter-clump medium model 
E (in dashed red line) to the observational data (in 
solid black line). Right: Comparison of the anisotropic 
model F (in dashed red line) to the observational data 
(in solid black line). Model parameters are given in 
Table 2. 



Fig. 19.— Even though the O VIII A18.97 A line is 
well fitted with an homogeneous model (Model A in 
dashed red line), the O VII triplet and N VII lines in 
the model are too weak. That indicates the need to use 
a second temperature in order to fit the whole spectrum 
of ^Puppis. Consequently, the hot gas filling factor of 
the first plasma temperature will be lower in order to 
take account of the contribution, in the O VIII inten- 
sity and line profile, of the second plasma temperature 
emission. 
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